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Electronic transport properties in a disordered quantum wire are very well described by the 
Dorokhov-Mello-Pereyra-Kumar (DMPK) equation, which describes the evolution of the transmis- 
^— J , sion eigenvalues as a function of the length of a multichannel conductor. However, the DMPK 

. equation is restricted to quasi one dimensional systems only. We derive a generalized DMPK equa- 

' tion for higher dimensions, containing dependence on the dimensionality through the properties of 

the transmission eigenvectors, by making certain statistical assumptions about the transfer matrix. 
^ . An earlier phenomenological generalization is obtained as a special case. 
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The Dorokhov-Mello-Pereyra-Kumar (DMPK) equation || has been enormously successful in describing the 
electron transport properties of a quasi one dimensional disordered conductor ||^ . The equation describes the evolution 
of the joint probability distribution of the transmission eigenvalues with increasing length of the system, and has been 
shown to be equivalent |^ to the description in terms of a non-linear sigma model obtained from the microscopic 
^ , tight binding Anderson Hamiltonian. The advantage of the DMPK approach over the non-linear sigma model is that 

a one can consider the full distribution of transport quantities rather than the mean and the variance alone. Recent 
analytical as well as numerical results show that the distribution of conductances has many surprises, including very 
sharp features at (dimensionless) conductance <? = 1 which could not be anticipated from studies of the moments of 
the distribution, and which should have important consequences for the Anderson transition. One major disadvantage 
of the DMPK equation, however, is that it does not contain information about the spatial structure of the sample in 
' directions perpendicular to the direction of the current flow, limiting its applicability to quasi one dimension (QID) 
, only. Since at present there is no other analytic approach available to study the full distribution of transport properties, 
Q ■ it is clear that a generalization of the DMPK equation valid in higher dimensions is of fundamental importance. A 
^ O ^ , phenomenological generalization, with an ad hoc constraint to conserve probability, was recently proposed which 
seems to agree with numerical results in systems beyond QID in some restricted regimes. In the present paper 
' we derive a further generalization that contains dependence on the dimensionality through the properties of the 
^ transmission eigenvectors, and contains the earlier model as a special case. No additional constraint is needed 
*n ' to conserve probability in the present approach. Moreover, the approach reproduces the expression for Lyapunov 
. exponents in higher dimensions obtained in Known properties of these exponents provide useful constraints on 
the phenomenological parameters in the current model. 

In the transfer matrix approach, a conductor of length L is placed between two perfect leads of finite width. The 
■ scattering states at the Fermi energy define TV channels. The 2N x 2N transfer matrix M relates the flux amplitudes 
\ on the right of the system to that on the left jl^ . Flux conservation and time reversal symmetry (we consider the 
"Y^ • case of unbroken time reversal symmetry only) restricts the number of independent parameters of M to N{2N + 1) 
a I and can be written in general as |^ 

1; -^ii ::.){'W J^){i :-)-^"^-' 

O ■ where u,v are N x N unitary matrices, and A is a diagonal matrix, with positive elements X.i,i = 1,2, ...A^. An 
ensemble of random conductors of length L, all with the same macroscopic disorder characterized by the same 
mean free path I but different microscopic realizations of the randomness, is described by an ensemble of random 
transfer matrices M, whose differential probability depends parametrically on L and can be written as dpL^M) — 



5L^ L 



FIG. 1: A small wire with length 5L, "the building block", is attached to a long one with length L {SL <^ L). The total transfer 
matrix M is given by M = M" M' . 
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PL{M)dfi{M) = pL(r, U, V)dfi{r)diJ,{U)dfi{V). Here dii{M) is the invariant Haar measure of the group, given in terms 

of the parameters in (1) by dfi{M) = J(A) dXi dfj,{u)dfi{v), where J(A) = Yl^^j \Xi ~ Xj\'^ with the 'symmetry 

exponent' /3 = 1 in our case, and d^{u) and d^iv) are the invariant measures of the unitary group U{N). 

We now take two conductors, with lengths L and SL, to make a conductor of length L + SL (Fig. 1). Denoting the 
corresponding transfer matrices by M", M' and M = M"M' with probability densities pl{M" = MM'^'^), psl{M') 
and Pl+sl{M), respectively, we have the relation [^j: 

Pl+sl{M)^ j PL{MM'-^)psL{M')dii{M'). (2) 

The restriction of the DMPK equation to QID arises from the "isotropy" approximation, that the distribution 
Pl{M) is independent of the matrices u and v. Several attempts have been made in order to relax the isotropy 
approximation |0, ^, We will avoid writing a set of coupled evolution equations for A, u and v by considering the 
marginal distribution 

PLix) = J pl{t, c/, v)d^l{u)d^l{v). (3) 

We first show that the U integral can be done exactly and then we make statistical assumptions about certain products 
of the matrix elements of V. 

Consider the combination H = MA/t = UT^W . At length L, we get H" = W'M'"^ = UTVM'~^{M'-^)W^TU^ = 
jj//p//2jj//t_ ^jjgjj follows that we can write U" ^ U ■ U"'{T,V,M'). Similarly, by considering the combination 
Q = MHI = V'^T^V, we get Q" = M"'^M" = {M'-^)'^V^T^VM'-^ = V'^T'^V", which implies that we can write 
r" = r"(r, V, M') and V" = V"{T, V, M'). Eq. (2) can then be rewritten as 

Pl+sl{U,T,V) = J pl{U ■U"'iT,V,M'),T"{T,V,M'),V"{T,V,M')) 

X psLiM')d^iiU')dii{r')d^iiV'). (4) 

We now integrate both sides of Eq. (4) over the invariant measure dfi{U). The left hand side defines a marginal 
distribution qL-^-SL(^,V). On the right side, since the measure is invariant, dpi{U) = diJL{U ■ U'") for U'" fixed. 
Therefore diJ.{U) — diJ.{U"), and integrating over diJ.{U") gives the marginal distribution (/^(r", V") with the following 
convolution equation: 

qL+5L{\v) = J qL{y',v")psL{X',u',v')dpi{M'), (5) 

where we have used Eq. (1) to introduce the matrix elements A, v, A', w', etc. 
Writing A" = A + SX and v" ~ v + Sv, Eq. (5) can be rewritten as 

qL+SL{>^,v) = {qL{X + SX{X,v),v + 5v{X,v)))g^ , (6) 

where {■■■) denotes an average over the ensemble of M' . In order to obtain SX and Sv within a perturbation theory, 
we search for a matrix constructed from M whose eigenvalues and eigenvectors are given by A and columns of v, 
respectively. Consider the matrix Q = M'^M = V'^T'^V. Flux conservation implies — E^QS^, where is the 
2N X 2N generalization of the Pauli matrix cr^. It then follows that the matrix X = [Q + Q^^ — 2/]/4, where / is 
the identity matrix, is block diagonal. It has been shown that V diagonalizes X leading to N doubly degenerate 
eigenvalues A. We can therefore obtain i5A and Sv by considering the change SX — X" — X due to change in M 
arising from M' . Writing SQ = Q" — Q, it is easy to see that SX is also block diagonal, with SXu = ■^SQu and 
SX22 = ^SQIi- Since both X and SX are block diagonal, the perturbation can be treated as acting separately on the 
two sub-blocks of X, and one can use ordinary, as opposed to degenerate, perturbation theory to obtain SX and Sv 
by considering one sub-block only. 

Denoting the perturbation by w = v{SXii)v^ = ■^v{SQii)v^ we get 

w = -X + vu' \X' + VTTXJv'v^ Xvv'^f VTTy + Vyv'*v*'' Xv*v'*'' y/y] u"'v'< 
- vu'[VTTX'v'v''y'X{l + X)v*v'*'< Vy + VX'v'*v*'< X{1 + A)WV(1 + A']u'1'u^ (7) 
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We expect A' oc 5L/1 « 1. Since w contains terms proportional to vA', we need to consider corrections to both 
A and v up to second order in w in order to keep terms up to O(A'). Standard perturbation theory [p^ gives 

6Xa = <5Ai'^ + <5Ai'^ and 6vl„ = Svl^n'^ + Svl^^'^ with 

a«=^.„; (8) 



^^t(i) ^ Y: Y^vL (9) 

and similarly for 6van ■ Note that v + 6v has to remain unitary, which imposes an additional constraint. 

The averages over the "building block" M' appearing in Eq. (6) involve averages over combinations of A', u' and 
v' which appear in w. Note that the building block is highly anisotropic; in the limit SL —^ 0, the transfer matrix 
M' —^ I. This condition will be implemented by assuming u'v' = I. In addition, instead of modelling the full 
Psl{X' , v'), we use the following averages over A/': 



a aa aj 




(10) 



ab I SL 



where k = SL/l. The first average is used in ||], the second one in y]; the latter incorporates the fact that the thin 
slice allows backward scattering without changes in the channel indices, and is highly anisotropic. Eqs. (10) and (11), 
with the conditions u'v' — I and k <C 1, define the model of our building block. 



Using the above model and the expansion \/l + A' = 1 + A'/2 + 0(A'^), we average SX, Eq. (8), over M': 

f Ab + 2, 

Aa — Afc 



{SXaSXb)sL = Sab2KXa{l + Xa) ^ \Vaaf ■ (13) 

a 

Similar calculations can be done to obtain the corresponding averages for 6v. 

We should now, in principle, expand Eq. (6) in a Taylor series around both A and v and evaluate the average over 
M' . However, for weak disorder, it is known that the isotropy approximation is very good, which means dqL /dvab — > 0. 
In the strong disorder limit, the eigenvectors are localized, so that in the the nth row Vna has only one element equal to 
unity, all other elements being zero. Since different rows are orthogonal to each other, we can see from (9) that in this 
limit 6v because of the restricted sum m ^ n. So in both these limits, the product (dqL/dvab) < Svab >sl—* 0. 
To a good approximation, in these limits, we can therefore treat v at the macroscopic L as a parameter that depends 
on disorder but for a given strength of disorder does not change any further with SL. In general, a conductor has a 
fraction (depending on the disorder) of its channels closed ||l^, i.e. the corresponding eigenvectors are localized, while 
the others are open, the corresponding eigenvectors being extended. We expect the isotropy approximation to remain 
good for the open channels, and the closed channels to contribute Sv ^ since they are localized. In other words, the 
product (dqL/dvab) < Svab >sl can be assumed small (compared to the contribution from {dqL/dX) < SX >sl) for 
all channels for a wide range of intermediate disorder as well. Physically, since the eigenvalues depend on the length 
exponentially, SX{X, v) remains important for all lengths at any disorder. On the other hand, we expect that at a 
macroscopic length L, the eigenvectors already evolve to either metallic (isotropic) or insulating (localized) structures 
for a given macroscopic disorder, and any further change due to SL (as opposed to change in disorder) is likely 
to be negligible. We will therefore expand the ensemble averaged marginal probability density ^^(A, w) within the 
approximation 

{qUX + SX{X, v),v + Sv{X, v))),^ « {qUX + ^(A, v),v)),^ . (14) 
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We show below that this approximation retains the dominant eigenvector correlations (via S\{X,v)) needed to repro- 
duce the Lyapunov exponents in arbitrary dimensions as obtained in . The price we pay for not including Sv in our 
calculation is that we will not be able to evaluate the eigenvector correlations self consistently, but will have to use 
them as phenomenological parameters. 

Using (14), we expand (6) in a Taylor series about A, 

a ah 

where the dots include terms containing higher order derivatives of A. We choose k = <C 1, which allows us to 
truncate the Taylor series at the third term, neglecting terms of 0{k^) and higher. Using Eqs. (12) and (13), the 
reduced distribution PlW = J qLiX,v)dfi{v) can then be written as 



Pl+sl (A) «pi(A) + KE(l + 2Aa) 



, -sr^ Xa + Xb + 2XaXb f ST^ . ,21 \2'^1l{\v) 

+ ^Y.^Xail + Xa) j E|.,„|4^!^l^d^(«). (16) 

a a ^ 

In order to make further progress, we will now make a 'mean-field' approximation, where the products of four w's 
that appear inside the integrals in (16) are replaced by their mean values which can be taken out of the integrals. 
This is equivalent to the assumption that for a given disorder, fluctuations in such quantities are small compared to 
their averages. In the weak disorder regime, each matrix element is of order l/\/]V, differing mostly in their phases; 
once the phases cancel, the fluctuations are negligible for homogeneous disorder. In the strong disorder regime, each 
eigenvector has one element which is unity representing a localized site and the rest are zero, but different samples 
will have different localized sites leading to large sample to sample fluctuations for individual elements. However, 
it is expected that the fluctuations in the sum over the elements of any eigenvector for a given disorder will remain 
negligible. Therefore the mean-fleld approximation is reasonable in these two limits. As argued before, we will assume 
that the assumption remains valid in the intermediate region of disorder as well, based on the picture of open and 
closed channels. Within this approximation, and expanding the left hand side of Eq. (16) in powers of 5L/1, we get 

dpLW , .dpLiX) y-^ Xa + Xb + 2XaXb dpL{X) 

a {LI) ^ dXa ^ Xa-Xb dXa 

+ EA.(1 + A.)i^..%|^. (17) 



Here we have defined 



J2{Ko.\^M^)l = Y1 J kaa|'k6a|V(A,«)dAi(A)dM(w) =ifafc. (18) 



Since v is unitary, Kab satisfies the sum rule J2b ^a.b = 1- This allows us to rewrite Eq. (17), following Q, as 



dpL{X) _ 1 V- J_ 



- dp 

d{L/l) J ^ dXa L^"^^ ^ ^"^^"°'^9A:, 
with 



(19) 



J=niAa-Afcr-, Jab=^. (20) 

a<b 
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Equation (19), with the definition (20), is our generahzation of the DMPK equation. Note that in the DMPK 
equation, time reversal symmetry fixes the exponent in the Jacobian J to be /3 = 1. In our case as the eigenvectors v 
are integrated over, the coupUng between the eigenvalues and eigenvectors adds an effective matrix exponent to the 
existing symmetry exponent, resulting in an effective J in eq. (19). 

Under isotropy condition Kab = ^^^i , we recover the DMPK equation {'jab = 1)- If we choose Kab — 
and Kaa ~ 'n+i^ we obtain the generalization of [7], where an extra condition was needed between /ii and /i2 in 
order to satisfy the conservation of probability. In our current framework, that condition is identical to the sum rule 

Eb^ab = l. 

As a check of our model, we evaluate the Liapunov exponents i^a = 23r ^^^^ X^) '^^ ^^'^ limit Xi ^ X2 ^ ■ ■ ■ ^ 
Xn S> 1 and compare with Expanding in powers of dXa/Xa, averaging over v and using the results for < 6X >, 
Eqs. (12-13), we obtain 

2l^a^^2^i + 2Kj2 (^It'aal'kban -K/^l^aalM • (21) 

bi^a) °^\q / L \ a I L 

Separating the sum over h into a term less than a and another greater than a and using the unitarity of w, Ua can be 
rewritten in the form 

1 ^ 

l^a « ^Kaa + Yl ^'''>, (22) 

b=a+l 

which coincides with that of Note that the isotropy approximation gives the smallest Liapunov exponent to be 
'^N = ivTT' ^^^^^ ^^^^ ^^"^ correct form in the metallic regime (Ai <C 1), the approximation fails to reproduce the 
expected behavior in the insulating regime, ~ 0(1) independent of N. Our approach incorporates the necessary 
eigenvector correlations to describe the Liapunov exponents at weak as well as strong disorder limits. Known properties 
of these exponents should provide a guide for constructing a model of K. 

It is important to note that while Eq. (19) is of the same form as the DMPK equation, the presence of the matrix 
K might not allow a solution of (19) in the same generic form as that of the DMPK equation ||^. 

In summary, the DMPK equation for the distribution (A) was obtained by adding a thin slice of length 5L to a 
conductor of macroscopic length L. The assumption that pl depends only on A and not on (u, v) limited the equation 
to QID only. We consider the marginal distribution where the u and v are integrated over. The u integral is done 
exactly. We then assume that while the changes in A due to the added slice depend crucially on the eigenvectors 
V, the eigenvectors themselves do not change much with length, i.e. they remain either metallic or insulating as 
already determined at length L. This implies that our equation is valid only beyond the relaxation length of the 
parameters Kab- We also assume that the eigenvector correlations \vaa\'^\vba\'^ have sharp distributions. While 
these assumptions were made to obtain the simplest generalization that captures the essentials of dimensionality 
dependence and need to be verified independently (e.g. numerically), we show that our approach keeps the dominant 
eigenvector correlations that reproduces the Liapunov exponents in higher dimensions at both weak and strong 
disorder. The fact that the parameters Kab may incorporate proper dimensionality dependence has already been 
shown in Finally, equation (19) reduces to the DMPK equation as well as to an earlier generalization in 
appropriate limits. 
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